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Abstract 



We investigate the temperature dependence of the instanton contents of gluon 

fields, using unquenched lattice QCD and the cooling method. The instanton 
L-j | size parameter deduced from the correlation function decreases from 0.44fm 

^^ ■ below the phase-transition temperature T c (~ 150MeV) to 0.33fm at 1.3 T c . 

The instanton charge distribution is Poissonian above T c , but it deviates from 
^ . the convoluted Poisson at low temperature. The topological susceptibility 

decreases rapidly below T c , showing the apparent restoration of the U(1)a 

symmetry already at T w T c . 
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Among the four global symmetries associated with quantum chromodynamics (QCD), 
two chiral symmetries, SU(Nf)A and U(1)a, must be broken as the observed hadron spec- 
trum suggests. At high temperature, it is expected that the two symmetries are restored. 
The question as to how they are restored is an important issue, especially in relation to 
the QCD deconfinement phase transition (or maybe a crossover), and is recently attracting 
much attention [!]-§]. The deconfinement phase transition is thought to have occurred in 
the evolution of the early universe and is expected to generate a new phase of matter, the 
quark-gluon plasma (QGP), which is under investigation in relativistic heavy-ion collisions. 

Instantons, tunneling between topologically different gauge vacua, are believed to play 
crucial roles in the breaking and restoration of the chiral symmetries, and perhaps also in the 
deconfinement phase transition |2[]: The SU(Nf)A symmetry manifests itself in the Goldstone 
mode through chiral condensates generated dynamically by instantons. A similar mechanism 
may be at work in the formation of hadronic bound states, and the SU(Nf)A restoration 
and the deconfinement are thus closely related to each other. The U(1)a symmetry is 
spontaneously broken by the chiral anomaly driven by instantons. At high temperatures, 
the restoration of the symmetries takes place because instanton amplitudes are suppressed 
through Debye screening, which induces the gluons to be (electrically) massive. 

How realistic is this instanton mechanism, and what is the precise dynamics of the in- 
stanton participation? The best means to answer the questions is the lattice gauge theory 
||, which provides a systematic, nonperturbative framework to study the equilibrium prop- 
erties of QCD at both zero and finite temperatures. Whereas direct comparison between 
lattice results and experiments is not always possible, lattice QCD often motivates and con- 
strains phenomenological models, which connect the QCD with experiments. Here, the use 
of unquenched QCD must be essential, since the instanton mechanism is expected to be 
affected by the presence of the dynamical quarks. 

In this letter, we report unquenched lattice QCD results on the temperature depen- 
dence of the instanton contents of gluons. For the first time, we find curious temperature 
dependence suggesting that instantons are seemingly correlated below the SU{Nf)A phase 
transition, T c « 150MeV. The nature and the details of the correlation are yet unclear, 
however, and so is the question of whether it corresponds to any of the dynamical models 
@. We also find that the U(1)a symmetry appears to be largely already restored at T c . 
Note that because of the cooling method we used, we observed neither small instantons nor 
instanton molecules. 

The calculation is carried out for 16 3 x N t lattices with N t = 4, 6, 8, 10, 12, 14, and 16 
for m q a = 0.0125. The temperature is varied by changing ]V t , with (3 = 5.54 fixed. The 
susceptibility and the topological-charge correlation are observed to be stable with respect to 
the number of cooling steps. We find that the topological susceptibility decreases slowly at 
low temperature as the temperature increases, but a rapid suppression sets in around T = T c , 
so that, compared to the value at T = 0, it is already down by an order of magnitude at T c 
and becomes only about 1% at 1.3T C . Accompanying this rapid change in Xt is a transition of 
the instanton size from p « 0.44fm at low temperature to 0.33fm at 1.3T C . The distribution 
of instanton charges is Poissonian above T c , but it deviates from the convoluted Poisson at 
low temperatures. 

The instanton content of the gauge fields is monitored by the the topological-charge 
density, which can be measured on the lattice with 



Q{x n ) = -—-^e a/3lS Re Tr [U af3 (x n )U lS (x n )] , (1) 

where U a p is the product of the link variables around a plaquette in the a — [3 plane. A direct 
calculation of Q from gauge configurations is plagued with poor statistics and with lattice 
artifacts associated with discretization. Teper et al. have shown that the topological 
charge can be reliably extracted using the cooling method, in which one applies a certain 
number of smoothing steps to each configuration in order to minimize the action locally. The 
ultraviolet fluctuations are quickly suppressed by these cooling steps, and only the relatively 
stable instantons are left in the configuration, making it easy to extract the topological 
charge. Details of the cooling method are discussed in |6|.[7|] . 

The topological susceptibility is then given as the fluctuations of the topological charge: 





* s ™o< £GM > , (2) 



where (...) indicates configuration averaging, N t and N s are the number of sites in the 
temporal and spatial direction, respectively, and a is the lattice spacing. 

In order to examine the details of the topological-charge distribution, we also calculate 
the topological charge- density correlation function [7|, 

C Q (x) = (J2Q(y)Q(x + y)}/(j:Q 2 (y)} . (3) 

y y 

One can compare this correlation function to a convolution of an isolated instanton 
topological-charge density 
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Q^) = -^\-r-r,\ • (4) 



Zero-temperature calculations using the cooling method show that after about 50 cooling 
steps, the gauge fields are dominated by large, isolated instantons, whose profiles agree well 
with Eq. (§. 

To generate the unquenched gauge configurations, we adopted and modified a code writ- 
ten by the MILC collaboration Hj, using its option of dynamical Kogut-Susskind fermions. 
Forty configurations are generated for each of N t = 12, 14, 16, and 100 configurations for 
each of N t = 4, 6, 8, 10. In each case, a molecular-dynamics time step of dt = 0.02 is 
used, and the first 300 time units (43 for N t = 12, 14, 16) from a hot start are used as 
thermalization. The configurations were separated by 3 time units for N t = 12, 14, 16; 5 
time units for N t = 4, 6, 10; and 10 time units for ]V t = 8. 

These configurations are then cooled, using the standard procedure as used in 0, for 200 
steps, saving every 20 steps. We have checked to see that both the topological susceptibility 
and the topological-charge correlation function settle to stable values by 100 cooling steps. 
The results reported here are based on those extracted after 120 cooling steps. 

The temperature of a gauge configuration is given by T = l/(N t a) and can be varied 
by changing either N t , or a ( = the lattice spacing), a is a function of the lattice inverse 
coupling j3, and can be calculated by the perturbation theory, but only for a narrow range 
of /3 (and hence to T) in the asymptotic scaling regime. In order to obtain the temperature 



dependence reliably, we choose to vary the temperature by varying N t , keeping a and j3 
fixed. In this way, the uncertainty in a enters only as an overall constant factor, which little 
affects the functional relation of the Xt and T. The price we pay is that the calculation 
can be done only for a limited number of temperatures corresponding to the discrete values 
of Nt- This is especially problematic at high temperatures. The phase transition with the 
parameters we use occurs at N t = 8 and a temperature of T c = 150(9) MeV [(J, giving 
a m 0.17fm. 

Figure 1 shows the topological susceptibility as a function of the temperature. The solid 
curve in the figure is the PCAC suppression that is due to soft pion gas [ID], which is believed 
to be relevant for low temperature T «T C : 



Xt(T) = X t(T = 0)(l + cT 2 /F;- 



(5) 



Equation |5| fits our data very well up to T = T c , but with c = —0.36, a value that is well 
outside the allowed range based on the PCAC argument — 1/6 < c < 1/6 JL(|. At higher 
temperatures, the suppression of the susceptibility follows a pattern different from the PCAC 
prediction. 

Also shown in Fig. 1 is the perturbative Debye-screening prediction (for N c = 3, Nf = 2) 



11 : 



Xt (T) = Xt(T = 0) (l + A 2 /3) 5 exp -8A 2 /3 - 14a (l + 7 A 



(6) 



where A = npT, a = 0.01289764, and 7 = 0.15858. The dot-dashed and dashed curves 
correspond to p = 0.44 and 0.33 fm, respectively. Equation |5] is not expected to be valid 
below T c , but should become relevant at high temperatures. We observe that the curve 
for p = 0.44 fm fits the four lowest temperature points quite well, giving an alternative 
suppression model to PCAC at low temperatures, and that the curve for p = 0.33 fm 
fits the high temperature points quite well. This observation suggests that the Debye- 
screening mechanism may be already operative for T < T c , though the dominant instanton 
size parameter changes suddenly across T c . Unfortunately, our N t = 4 (T = 2T C ) result is 
rather uncertain because of severe finite-size effects and because not one single instanton is 
found in our 100 configurations. Also, our method of temperature variation by varying Nt 
provides only sparse data in the high-temperature regime. The temperature dependence of 
the Debye-screening predictions must be explicitly verified by further calculations. 

The general features of the temperature dependence of the Xt are the same in the 
quenched calculation [§] and in the present unquenched calculation, but detailed features 
differ significantly. The unquenched Xt decreases more slowly at low temperature, but goes 
down much more rapidly around T c . The Xt is suppressed by almost 90 % at T = T c in the 
unquenched case, while it is only 50 % in the quenched case. That is, the suppression of the 
instanton amplitudes is enhanced by the dynamical quarks. The rapid decrease in Xt con- 
tinues above T c , so that only about 1% remains at 1.3T C in the unquenched, and about 10% 
in the quenched calculation. This temperature dependence of Xt suggests that the window 
of temperatures above T c in which £7a(1) is broken is narrowed by the presence of dynamical 
quarks, and that Ua(X) is practically restored at about T c . One must take caution, how- 
ever, in applying lattice results on the U(1)a symmetry, because the symmetry is explicitly 



broken in the conventional lattice and the process of approaching the continuum limit in- 
volves uncontrolled subtlety ||. On the basis of the present calculation with m q a = 0.0125, 
we draw no conclusion on the relationship between T c and the U(1)a symmetry breaking 



temperature, and thus on whether the phase transition is of the first or second order [|T2| , [3~3 

Figure 2 shows the temperature dependence of the topological charge correlation func- 
tions in comparison to the profile of an isolated instanton, Eq.((|). In the quenched calcula- 
tions, Cq(t) is practically independent of the temperature, following remarkably closely to 
the profile, though hints of a transition to a smaller size parameter have been observed above 
T c j8|. The present unquenched calculation confirms the same trend below T c , but it clearly 
shows a rapid transition to a smaller size around T c . For clarification, we show in Fig. 2 the 
temperature dependence only around T c : the size parameter, p, is approximately 0.44 fm at 
0.75T C , shrinks rapidly around T c , and becomes about 0.33fm at 1.3T C . It is interesting to 
note that the low-temperature value of p agrees fairly well with that used in the instanton 
liquid model, while the high-temperature value agrees with the dilute-gas model and also 
with the quenched value [0. 

The distribution of instanton charges in our configuration samples is shown as histograms 
for each N t in Fig. 3. Whereas in the quenched case, these distributions agree with the con- 
volued Poisson distribution || , the unquenched results show significant deviations from this 
Poissonian, with curious suppression of the distribution around the zero topological charge 
(Q = 0) for all Nt below the phase transition. The dynamical quarks have apparently in- 
duced non-trivial correlations among the instantons at low temperatures. Above the phase 
transition, the distributions become Poissonian, suggesting a return to the dilute-gas ensem- 
ble. 

The instantons we have discussed so far in this letter are all of the topological charge 
of the unit magnitude. We note that we have observed occassionally the instantons of the 
topological charge of the magnitude two, which are rarely treated in any dynamical model. 

In conclusion, our investigation of the temperature dependence of the instanton con- 
tents of unquenched gluon fields yields the following results: The topological susceptibility 
decreases rapidly around the phase-transition temperature T c « 150MeV; at T = T c , it 
is suppressed by one order of magnitude compared to its zero-temperature value, and at 
T = 1.3T C , it is reduced by another order of magnitude. The temperature dependence of the 
topological susceptibility is consistent with a PCAC prediction for temperatures up to T c , 
though with a large coefficient c. It also agrees with a Debye-screening model, but with a 
sudden shrinking of the instanton size parameter at around T c . The topological charge cor- 
relation function also changes rapidly around T c , which can be characterized by a transition 
of the dominant instanton size parameter from 0.44fm at low temperatures to 0.33fm at 1.3 
T c , The distribution of the instanton charges suggests the existence of nontrivial correlations 
among the instantons at low temperatures, which disappear above T c . Our results indicate 
that the Ua(1) symmetry is apparently restored largely at T « T c . 
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FIGURES 

FIG. 1. Topological susceptibility as a function of temperature. The data points corre- 
spond to our unquenched-QCD-plus-cooling results with N t = 6,8, 10, 12, 14, 16, and a = 0.17fm 
is assumed. The transition temperature T c is approximately 150 MeV for our calculations with 
P = 5.54, ma = 0.0125, iVa; = 16. The dashed curve is the prediction of a PCAC model, Eq. (|5|), 
with c = —0.36; and the solid and dot-dashed curves of a Debye-screening model, Eq. (||), with the 
instanton size parameter, p = 0.44 and 0.33 fm, respectively. 

FIG. 2. Topological charge correlation functions across the phase-transion temperature 
(T c = 150MeV) from the same calculations as in Fig. 1. 

FIG. 3. Histograms showing the distribution of topological charge in the configurations for 
Nt = 6, 8, 10, 12, 14 and 16. Q is the magnitude of the number of topological charges in a configu- 
ration, and N(Q) is the frequency. The dots in the figure are the convoluted Poisson distributions 
fitted to the histograms. 
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